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Abstract 

For the first time, we obtain the complete and exact analytical solution of 
the Lawrence-Doniach model for layered superconductors in external paral- 
lel magnetic fields. By solving a nontrivial mathematical problem of exact 
minimization of the free-energy functional, we derive a closed, self-consistent 
system of mean-held equations involving only two variables. Exact solutions 
to these equations prove simultaneous penetration of Josephson vortices into 
all the barriers, yield a completely new expression for the lower critical held, 
refute the concept of a triangular Josephson vortex lattice and clarify the 
physics of Fraunhofer oscillations of the total critical Josephson current. 

PACS numbers: 74.80.Dm, 74.20.De, 74.50.+r 

In this paper, we obtain for the first time the complete and exact analytical solution of 
the popular phenomenological Lawrence-Doniach QXj (LD) model for layered superconductors 
in external parallel magnetic fields. 

At present, there is a universal belief that the LD model can adequately describe high-T c 
superconductors exhibiting the intrinsic Josephson effect. [0 Surprisingly, despite a large 
number of publications on the LD model over the recent years, it has not been realized yet 
that in the presence of a parallel magnetic field the LD free-energy functional provides a 
rare example of exactly solvable models in theoretical physics. (As an exception, there is a 
particular exact solution of Theodorakis, || valid in a restricted field range.) Moreover, even 
a closed, self-consistent system of mean-field equations for the LD functional has not been 
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obtained up until now, which in some cases leads to spurious results. Thus, calculations of 
the lower critical field are based on an arbitrary assumption of single Josephson vortex 
penetration and a continuum approximation, incompatible with the discrete nature of the 
LD model. (This treatment was criticized in Ref. 3.) The claim |7]] that the Fraunhofer 
pattern of the critical Josephson current occurs in the absence of Josephson vortices is at 
odds with the well-known situation || in a single junction. Furthermore, a hypothesis of a 
triangular lattice of Josephson vortices stands in direct contradiction to the exact solution of 
Ref. 3. Here we show that the origin of these inconsistencies is in an incorrect mathematical 
approach to the minimization of the LD functional, neglecting principal aspects of any gauge 
theory. [|10j Based on exact variational methods, we derive a remarkably simple, closed, self- 
consistent set of mean-field equations involving only two variables. As these equations 
turn out to be a special limiting case of a recently developed microscopic theory, |Ll| we 



concentrate here mostly on the problem of exact minimization of the LD functional and 
provide a brief summary of the main new physical results at the end of the paper. 

We begin by reminding basic features of the LD model. [fL,|12" In this model, the tem- 
perature T is assumed to be close to the "intrinsic" critical temperature T c0 of individual 
layers: 
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The superconducting (S) layers are assumed to have negligible thickness compared to the 
"intrinsic" coherence length C(T) oc t -1 / 2 , the penetration depth A(T) oc t -1 / 2 , and the 
layering period p. Taking the layering axis to be x, choosing the direction of the external 
magnetic field H to be z [H = (0,0, H)} and setting H = c = 1, we can write the LD 
free-energy functional as 
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Here A = (A^A^O) is the vector potential, continuous at the S-layers: A(np — 0,y) = 
A{np + 0,y) = A(np,y); W z is the length of the system in the z direction; f n (y) [0 < 
fniy) < 1] an d 0n(y) are, respectively, the reduced modulus and the phase of the pair 
potential A n (y) in the nth superconducting layer: 



A n (y) = A(T)f n (y)exv<f) n (y), 
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with A(T) being the "intrinsic" gap [A(T) oc r 1 ' 2 ]; H C (T) is the thermodynamic critical field; 
r(T) = 2a p hT~ 1 is a dimensionless phenomenological parameter of the Josephson interlayer 
coupling (0 < a p h 1). The local magnetic field h = (0, 0, h) obeys the relation 

dAy(x,y) dA x (x,y) 
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Our task now is to establish a closed, complete, self-consistent system of mean-field 
equations of the theory, which is mathematically equivalent to the minimization of @ with 
respect to /„, n , and A. First, we want to point out a common mistake |13| in the approach 
to this problem: It has not been realized in the literature that variations with respect to 
<p n and A are not independent and do not yield a complete set of equations. Indeed, as the 
functional (0) is invariant under the gauge transformations 



4>n{y) -> <t>n(y) + 2eA(np, y), A^x, y) -> A^x, y) + <9jA(x, y), i = x,y, 
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where X(x, y) is an arbitrary smooth function of x, y in the whole region (— oo < x < +oo) x 
(L y i < y < L y2 ), variational derivatives with respect to <fr n , and A x> A y are related by the 
fundamental identities 



2e- 



5Q 



LD 



d 5Q 



LD 



+ 



5Q 



LD 



5Q 



LD 



(6) 



54> n (y) dy5A y (np,y) 5A x (np + 0,y) 5A x {np -0,y)' 

Being a consequence of Noether's second theorem, such identities are typical of any gauge 
theory. |10| They imply that the number of independent Euler-Lagrange equations is less 



than the number of variables, and complementary relations should be imposed to elimi- 
nate irrelevant degrees of freedom and close the system mathematically Whereas in bulk 
superconductors and single junctions the elimination of unphysical degrees of freedom is 
accomplished by fixing the gauge, in periodic weakly coupled structures this problem has 
additional implications. fTTJ Namely, in the presence of the Josephson interlayer coupling 
the quantities 3> nn _i are not independent but subject to a set of constraint relations. Un- 
fortunately, this fundamental feature was not noticed in any previous publications on the 
LD model. 

Varying with respect to A x , A y in the regions (n — l)p < x < np under the assumption 
SA x (x, Lyi) = 5A x (x,L y2 ) = yields 
dh{x ? y\ 
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where j n ,n-i{y) is the density of the Josephson current between the (n — l)th and the 
nth layers, j = r(T)p/167re( 2 (T)\ 2 (T). Minimization with respect to A y (np,y) leads to 
boundary conditions at the S-layers: 
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Equations (0)-© should be complemented by boundary conditions at the outer interfaces 
y = L y \Ly2. As we do not consider here externally applied currents in the y direction, the 
first set of boundary conditions follows from the requirement [jny] y=Lyl Ly2 = 0: 
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Applied to Eqs. ([]), these boundary conditions show that the local magnetic field at the 
outer interfaces is independent of the coordinate x: h(x,L y i) = h(L y i), h(x,L y2 ) = /i(L y2 )- 
The boundary conditions imposed on h should be compatible with Ampere's law h(L y2 ) — 
h(Lyi) = 4nl obtained by integration of Eqs. (0) over y, where 

Ly2 Ly2 
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is the total current in the x direction. 

Differentiating ([]) with respect to y and employing (0), we arrive at the current- 
continuity equations for the S-layers: 
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Adding Eqs. (|T2"D , integrating and using boundary conditions (0), we get the first integral 
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This equation has mathematical form of a constraint relation and states that the total 
current in the y direction is equal to zero. 

The Euler-Lagrange equations for (p n do not yield anything new and only reproduce Eqs. 
(|12|), as expected by virtue of Noether's identities (|6|). To obtain complementary constraint 
relations, closing the system of the Euler-Lagrange equations and minimizing the free energy, 
we must modify the variational procedure. 

Noting that the kinetic energy of the intralayer currents in (|2[) can be minimized inde- 
pendently of the Josephson term, we impose additional constraints 
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compatible with boundary conditions (0) and constraint relation (|13|). The requirement 
of compatibility with the current-conservation law ([12]) automatically yields another set of 
constraints 



^n.n— 1 (y) = fn+i(y) sin 

(y)- 
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The physical meaning of Eqs. (|14D and (|T5|) that provide the sought necessary conditions for 
the true minimum of the free-energy functional (0) is obvious. Constraints ([14]) minimize 
the kinetic energy of the intralayer currents (it proves to be identically equal to zero) and 
assure the continuity of the local magnetic field at the S-layers: h(np + 0, y) = h(np — 0,y). 
[See Eq. (§)]. These constraints appear already in the case of decoupled S-layers. On the 
other hand, constraints ( |i~5l) are uniquely imposed by the Josephson interlayer coupling. 



Their function is to make the Josephson energy stationary with respect to variations of < 
and to assure the continuity of the Josephson current at the S-layers. 

As no other conditions are imposed on the variables, we can satisfy fllED by choosing 



f n {y) = fn-tiy) = f(y), ^n+lM = ®n,n-l(y) = <%)• 
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The establishment of constraints (|14] ) - (|T6|) , minimizing the free energy and closing the set of 
mean-field equations, is a key result of this paper. For example, these constraints automati- 
cally rule out any possibility of previously proposed [Q single Josephson vortex penetration 
and the hypothesis of a triangular Josephson vortex lattice. It should be noted, how- 
ever, that both the exact solution of Theodorakis for the dense vortex state and early 
calculations of the upper critical field are fully compatible with relations (|T^)-(|TB|). 

The remaining unphysical degree of freedom, related to the gauge invariance, is elimi- 
nated by fixing the gauge: 



A x (x,y) = 0, A y (x,y) = A(x,y). 



(17) 



[Note that dA/dx and d 2 A/dxdy are continuous at the S-layers by virtue of (|9|), (0), and 
(0); (P~5D-] The second set of relations (|IB|) now yields 4> n (y) = n<p(y) + i](y), where 4>{y) is 
the coherent phase difference (the same at all the barriers), and r](y) is an arbitrary function 
of y that can be set equal to zero without any loss of generality. 

From (^|), using the continuity conditions for A, dA/dx and constraints flU]), we obtain 
A(x, y) = 2^^^, while the functional (0) becomes 
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where W x = L x2 — L x \. Minimizing (|TE|) with respect to f(y) [with arbitrary 5f(L y i), 
5f(L y 2)} and <p(y) [with 84>{L y i) = <50(L y2 ) = 0], we arrive at the desired closed, self- 
consistent set of mean-field equations 



A n (y) = Af(y) exp [in(f)(y)} , 
6 



(19) 



f(y) + ( 2 (T)^M- - f(v) - r{T) [1 - cos0(y)] f(y) = 0, (20) 
f(^) = f(^) = 0, (2D 

^ = ^T sm ^ ) ' (22) 

A J = (87rej p)" 1/2 , (23) 

= jn,n-i(y) = Jo/ 2 (2/)sin0(y) = i-^M (25) 
that should be complemented by appropriate boundary conditions on h(y) (see above) with 

Ly 2 

I = J dyj(y) , where j(y) is the density of the Josephson current. 

Remarkably, the coherent phase difference <j) (the same for all the barriers) obeys only one 
nonlinear second-order differential equation ( p2]) with only one length scale, the Josephson 
penetration depth Aj [Eq. (|2"3|)1, as in the case of the Ferrell-Prange equation for a single 
junction. [Mathematically, equation (|22] ) is a solvability condition for the Maxwell equa- 



tions.] Due to the factor f 2 , equation ( p2]) is coupled to nonlinear second-order differential 
equation ( pH|) describing the spatial dependence of the superconducting order parameter / 
(the same for all the S-layers). Equations (pi]) constitute boundary conditions for (p0|). The 
Maxwell equations (|24[) , (f25|), combined together, yield Eq. (|22|) , as they should by virtue 
of self-consistency. 

It is instructive to compare the above equations with those of previous publications, based 
on an incomplete minimization procedure. Thus, for $ n+ i >n ($/) one introduces |||7|,[13| an 
infinite non-self-consistent set of the so-called " difference-differential" equations, containing 
two length scales. By virtue of the constraint relations (fHf ), in the gauge fll7D this set reduces 



to only one equation (p2|) with f n (y) = 1, while the second length scale, \/2£(T)/ Jr(T), 
related to unphysical degrees of freedom, disappears from the theory. 
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On the other hand, equations (p!9|)-(|25|) are only a limiting case of the true microscopic 
equations, [|TTJ if one identifies r(T) with the microscopic parameter a£ 2 (T)/a£ and sets 
a/p = 0, where £ is the BCS coherence length, a is the S-layer thickness [£ <C a "C 



C(T),A(T)], and ol — 7/V3') / dttD(t) <§C 1 [D(t) is the tunneling probability of the barrier 
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between two successive S- layers]. 

Equations (|20|)-(|25|) admit exact analytical solutions for all physical situations of interest. 
Aside from the region near the second-order phase transition to the normal state (because of 
the unphysical assumption of negligible S-layer thickness, the LD model does not adequately 
describe this regime ||), these solutions stand in a one-to-one correspondence with those of 
the microscopic theory. |TT| For this reason, we only briefly summarize the main physical 
results here, accentuating differences between the exact solutions and previous non-self- 
consistent calculations. 

The local magnetic field is independent of the coordinate in the layering direction. [See 
Eq. fl21|).] The Meissner phase in semi-infinite (along the layers) samples persists up to the 
superheating field H s = (epAj) -1 . Contrary to previous suggestions, J3J Josephson vortices 
penetrate all the barriers simultaneously and coherently, forming peculiar structures that 
we term JTT| "vortex planes". The existence of a single vortex plane in an infinite (along the 
layers) sample becomes energetically favorable at the lower critical field i? c ioo = 2(7repAj) _1 . 
(Previous calculations [|]-[(| of iJ c ioo, based on an invalid assumption of single Josephson 
vortex penetration and an anisotropic continuum approximation, are incorrect.) In the 
fields H cloo <C H <C [ep(iT)]^ 1 , with r(T) <C 1, equations (|20| ) - (|25|) reproduce the vortex- 
state solution of Theodorakis. || (The triangular vortex lattice PJT5[] is not allowed by 
the exact equations.) The magnetization in the vortex state exhibits distinctive oscillatory 
behavior and jumps as a result of vortex-plane penetration. For a certain field range, our 
calculations yield a small paramagnetic effect. In contrast to previous assertions, J7j the 
Fraunhofer pattern for the total critical Josephson current in layered superconductors with 
W <C A j (W = Ly 2 — L y i) occurs due to successive penetration of the vortex planes and their 
pinning by the edges of the sample. The first zero of the Fraunhofer pattern corresponds 
to the lower critical field H c \w = n/epW of a finite sample. Finally, for the upper critical 
field H c2oo in an infinite layered superconductor, equations (|20D -([25|) yield the well-known 



S 



results, 



as expected. We conclude by observing that the established relation to the 



microscopic theory |TT[ casts light on the exact domain of validity of the LD model. 
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